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We give the governing equations for multiple scalar fields in a flat Friedmann-Robertson- Walker 
(FRW) background spacetime on all scales, allowing for metric and field perturbations up to second 
order. We then derive the Klein-Gordon equation at second order in closed form in terms of gauge- 
invariant perturbations of the fields in the uniform curvature gauge. We also give a simplified form 
of the Klein-Gordon equation using the slow-roll approximation. 

PACS numbers: 98.80.Cq JCAP03(2007)004, astro-ph/0610864v5 

I. INTRODUCTION 

In the cosmological standard model the universe undergoes at very early times a period of inflation, during which 
fluctuations in the scalar fields are stretched to super-horizon scales and later on source the Cosmic Microwave Back- 
ground (CMB) anisotropies and the large scale structure. In standard inflation models the scalar field(s) responsible 
for the accelerated expansion of the universe also provide(s) these fluctuations [1]. If we are interested in calculat- 
ing observational consequences of the scalar field dynamics, linear or first order cosmological perturbation theory is 
sufficient, if our main interest is in calculating the power spectrum of the fluctuations. However, if we are interested 
in higher order observables, such as the bispectrum, we also need cosmological perturbation theory beyond linear 
order. The study of higher order effects promises to give new insights into the physics of the early universe. Much 
progress has been made recently to provide the necessary theoretical tools to study inflation beyond the standard 
linear approximation. In this article we use second order perturbation theory to derive the Klein-Gordon equation 
governing the evolution of the scalar fields during and after inflation. 

At the moment there are mainly two approaches to study higher order effects and in particular non-gaussianity: 
One approach uses nonlinear theory and the separate universe approximation, either employing a gradient expansion 
[2, 3] or using the AN formalism [4-8] (see also [9] for early works). Both the gradient expansion and the AN 
formalism are so far only usable on large scales, that is scales much larger than the particle horizon, which is sufficient 
for many applications. The other approach uses second order cosmological perturbation theory following Bardeen 
[10-30]. The Bardeen approach is valid on all scales, it can be however more complex than the AN formalism. Both 
approaches give the same results on large scales. 

In this article we focus on deriving the Klein-Gordon equation at second order valid on all scales for multiple 
scalar fields minimally coupled to gravity in a flat Friedmann-Robertson- Walker (FRW) background spacetime. We 
take into account metric perturbations up to second order, considering scalar perturbations only. Throughout this 
article we shall work in the flat gauge. Having derived the field equations up to second order in the perturbations, 
we calculate the Klein-Gordon equation at second order in closed form for multiple fields. This allows to show that 
the Klein-Gordon equation also at this order does not contain terms that diverge in the large scale limit. These 
divergences might have appeared through the introduction of inverse gradients, but the calculations show that any 
inverse gradient present is balanced by one or more gradient terms. 

The paper is organised as follows: in the next section we give the governing equations, that is Einstein's field 
equations and the Klein-Gordon equation up to second order in the flat gauge. In Section III we give the Klein- 
Gordon equation at first and second order in closed form, that is solely in terms of the field fluctuations in the flat 
gauge, in real space. In Section IV we rewrite the second order Klein-Gordon equation in Fourier space. In Section V 
we employ the slow roll approximation to simplify the Klein-Gordon equation in Fourier space. We conclude in the 
final section. 

Throughout this paper we assume a spatially flat FRW background spacetime and use conformal time, r\. Derivatives 
with respect to conformal time are denoted by a dash. Greek indices, /i, v, A, run from 0, ... 3, while lower case Latin 
indices, i, j, k, run from 1, ... 3. Upper case Latin indices, /, J, K, denote different scalar fields. 
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II. GOVERNING EQUATIONS 

In this section we briefly review the governing equations of a system of multiple scalar fields. A more detailed 
exposition can be found for example in Ref. [29] . The covariant Einstein equations are given by 

G>„ = 8ttG , (2.1) 

where G^ u is the Einstein tensor, T^ v the total energy-momentum tensor, and G Newton's constant. Through the 
Bianchi identities, the field equations (2.1) give the local conservation of the total energy and momentum, 

V„T" , ' = 0, (2.2) 

where V M is the covariant derivative. The energy momentum tensor for N scalar fields minimally coupled to gravity 
is 

N 

I 

-9nvU(<pi,...,ip N ), (2.3) 



1 ad 
fK,ii^>K,v ~ -j9y,v9 <PK,a<PK,l3 



K=l 

where ipx is the Kth scalar field and U the scalar field potential and ipx.fj. = ^frr- 

Tensorial quantities are separated into a background value and perturbations according to 

T = T + 5Ti + i(5T 2 + . . . , (2.4) 

where the background part is a time dependent quantity only T = T (r/), whereas the perturbations depend on 
time and space coordinates x^ = [rj, x l ], that is ST n = 6T n (x^). The order of the perturbation is indicated by a 
subscript, i.e. <5T„ = 0(e n ), e here being the small parameter. 

The metric tensor up to second order, including only scalar perturbations, and without gauge constriction is 

5 oo = -a 2 (l + 2</> 1 +</> 2 ) , (2.5) 

5o* = a 2 (#i + \B^j , (2.6) 

9ij = a 2 [(1 - 2^i - V2) Sij + 2E hij + E 2iij ] , (2.7) 

where a = a(rj) is the scale factor, r] conformal time, Sij is the flat background metric, (f>i and 02 the lapse functions, 
and tpi and "02 the curvature perturbations at first and second order; B\ and B 2 and E\ and E 2 are scalar perturbations 
describing the shear. The contravariant form of the metric tensor is given in Appendix C 1. We shall employ the flat 
slicing and threading throughout this paper. In this gauge the spatial 3-hypersurfaces are flat and we have 

^ = ^ 2 = Ei = E 2 = . (2.8) 

Details on how to construct gauge-invariant variables at second order in general and in the flat gauge in particular 
can be found in Refs. [20] and [29]. We simply state the result here. 

The field fluctuation on uniform curvature hypersurfaces at first order, also known as Sasaki-Mukhanov variable 
[31, 32], is given for the Ith field as 

= Sipu + • (2-9) 

At 

The field fluctuation in the flat gauge, or Sasaki-Mukhanov variable, at second order for the Ith field is [20, 29] 1 



6 V2I = 5 V2I + ^ 2 + 



2 r 



2H(p' 0I + <p'oj - ij^foi 



+ 2|^Vi + ^i<W - 2^1,^ + X E) , (2.10) 



where X (tp, E) contains terms quadratic in gradients of the metric perturbations ipi and E\. Since we are working in 
the flat gauge throughout this paper, we shall drop the "tilde" to denote variables that have to be evaluated on this 
gauge. In particular, we denote the gauge invariant field fluctuations in the flat gauge simply by 5ipi and Sip 2 , hoping 
that this notation will make the equations more intuitive and transparent. Note however, that the field fluctuations 
in the flat gauge were denoted Qi and Q 2 at first and second order, respectively, in Ref. [29]. 



1 Rewritten using results from K. A. Malik and D. R. Matravers, arXiv:0804.3276 [astro-ph]. 
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A. Klein-Gordon equation 

We split the scalar fields ipi into a background and perturbations up to and including second order according to 
Eq. (2.4), 

tpiix") = ipoi(v) + Spuix") + • (2-11) 

The potential U = U(ipi) can be split similarly according to 

U{ m ) = U + SU! + \8U 2 , (2.12) 

where 

6U X = J2u VK 6<p 1K , (2.13) 

K 

5U 2 = J2 U 

K,L K 

and we use the shorthand U. VK = -§^- The energy-momentum tensor, Eq. (2.3), expanded up to second order in 
the perturbations for the metric tensor Eq. (2.5) is given in Appendix C2 (in the flat gauge). 

The evolution of the scalar fields is governed by the Klein-Gordon equation, which is the energy conservation 
equation, (2.2), for the energy momentum tensor defined in Eq. (2.3). Using the energy momentum tensor given in 
Appendix C 2 and the energy conservation equation (2.2) we get the Klein-Gordon equation in the multi field case at 
zeroth order for the Ith field 

^ I + 2H<p' 0I + a 2 U, Vl =0, (2.15) 

where H = at first order 

Sipu" + ZHStpu' + 2a 2 U Vl fa - V 2 5 VlI - if' 0I V 2 B 1 - ^\ + a 2 ^ U m 5<p 1K = , (2.16) 

K 

and we finally get at second order 

8if 2 i" + 2HSip 2 i' - V 2 (5^2/ + a 2 X! U ,<pi<pk<><P2K + a 2 U. ipiipKipL 5f 1K 5f 1 L + 2a 2 U. Vl (j)2 - ip' 0I (V 2 S 2 + <j>' 2 ) 

K K.L 

+ ^jBi,^* + 2 (2H<p' 0I + a 2 U. Vl ) B XM B Y k + 4^ (^a 2 ^ U, VIVK 6<p 1K - V 2 <Vi/^ + 4^j0i# 

- 28tp 1 i , (V 2 B 1 +<t/ 1 )-46tp 1 / ik B 1 * = 0, (2.17) 
where all the variables are in the flat gauge, as stressed in the previous section. 



B. The field equations 

In this section we give the field equations as given from Eq. (2.1) in the background, at first and at second order in 
the perturbations. Again we would like to stress that the perturbed equations are in the flat gauge throughout this 
section. 

At zeroth order we get the Friedmann equation from the — component of the Einstein equations 




(2.18) 



and from the i — j component we find 

(£) 2 - 2^ = SnG (j2 \p',k a'Uo) ■ (2.19) 
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1. First order 



Using the background equations given above, we find the — component of the Einstein equations at first order 
to be 



2a 2 U ^ + Vox W + a? Mi + T^V 2 ^ = 



K 



4ttG 



the — i part gives 



H<tn-4KG'£ i <f/ OK 5<p 1K = 0. 



K 



From the i — j component of the Einstein equation we get the trace free part 

B[ + 2HB 1 +<jn=0. 

Using Eq. (2.22) and the background field equations, we get the first order trace in its simplest form 



+ 4ttG 



a 2 (Wi + 2a 2 [/ o 1 - <Pok S( Pik' 



K 



= 0. 



(2.20) 



(2.21) 



(2.22) 



(2.23) 



2. Second order 

The — Einstein equation at second order, using the first and zeroth order — equations, Eqs. (2.18) and (2.20), 
is given by 



87TGa 2 U (02 + Bi, fc Bif) + HV 2 B 2 + 1 [b 1m B^ - (V 2 Bi) 2 ] - 2U^, k B k 

+4itGJ2 [4W + °? SU 2 + 4a 2 ^f/i<Ai + 6<p 1K ' 2 + 5<p 1K , k 8<pi K , k ] = .(2.24) 

K 

The — i Einstein equation is given by 

Ufa, - W.<h<h,i + mB 1M B£ + B 1M ^ k - V 2 B 1 ^ 1 . l - AttG^2 Wok <W,i + 2<W'<W,i] = . (2.25) 

K 

We then use the first order — i equation, Eq. (2.21), rewrite Eq. (2.25) and take the trace, which gives 
H (<j> 2 - 20? + B lik B^) - A-kG^^k^k 

K 

+V~ 2 (fa.kiB^ 1 - V 2 SiV 2 0i) -8^G^V- 2 (6<pi K 'V 2 6tpiK + 6tp' 1Kil 5<p 1Kt l ) =0, 

K 

where we introduce the inverse Laplacian, V~ 2 (V 2 )X = X. The i — j Einstein equation is given by 

|l67rGa 2 [/ (02 - 40 2 + B 1>k B^) + 2H<j>' 2 - 8H010; - 20i ifc ^* + AHB' l k B 1 k 

+V 2 (B 2 + 2HB 2 + <j> 2 ) - 20' 1 V 2 5i + B 1M B^ - (V 2 Bi) 2 + 32ttG {a 2 8U x + 2a 2 [/ o 1 ) 0i 



+8ttG 



X] (^ix,(^ix, ' ^ v'ok^k' - SipiK 2 ^ + a 2 SU 2 [Pj 

K J J 



(2.26) 



(2.27) 



(B' 2 + 2HB 2 + fa)*. + 20^ cf> hj + 2B^- (0' x + V 2 Bi) - 2B^ k B x k . - 167tG^<Vik ^(Wj = . 



K 



■5 



Taking now the divergence of Eq. (2.27) twice 2 we get 

8irGa 2 U Q (fa - H\ + B hk B k ) + Hfa 2 - ±Hfafa x ~ fa,kfa k + mB' 1>k B lt k + \ (b 1m B^ 1 - (V 2 Bi) 2 ) 



-fa 1 V 2 B 1 + 16ttG (a 2 J[/i + 2a 2 U fa) fa + 4ttG 



+V~ 2 {fa ! i fa, j +B i ; j (fa+V 2 B 1 )-B h i k B 1 k j -ZirGY J ^K, i 5yiK, 3 } 

K 



(2.28) 



where we made use of the relation [(/ 2 + f l9l ) <f J = V 2 (f 2 + f l9l ) = V 2 / 2 + /iV 2 ji + <?iV 2 h + ^./i The 
spatial trace is finally given by contracting Eq. (2.27) as 

24irGa 2 U (fa - H\ + B hk B k ) - 2fa M fa k - 2fa 1 V 2 B 1 - &HB k (2HB hk + fa tk ) 

+mfa 2 - UHfafa + V 2 (B' 2 + 2HB 2 + fa) + 1 [b^hB^ 1 - (V 2 Bi) 2 ] + iSirGfa (a 2 8U! + 2a 2 U fa) 

+4irG (3a 2 SU 2 - 3<p' 0K 5ip 2K ' - 3<W' 2 + <W,;<W, l ) = . (2.29) 

K 

In the above we made copious use of the background field equations without explicitly stating it. 



III. THE KLEIN GORDON EQUATION IN CLOSED FORM 

Using the field equations given in the previous section we can now rewrite the Klein Gordon equations at first and 
second order in the multiple field case. At first order the Klein-Gordon equation is, from Eq. (2.16), 



5<p u " + 2H5cp u ' - V 2 ^ 17 +a 2 ^ {u, VKVI + ^ (y>' iU, VK + <p' 0K U, Vl + ^o^o/^^o) } <W = ■ ( 3A ) 

Although there are notationally more compact forms of the Klein-Gordon equation at first order [33, 34], we find 
this form particularly easy to use. All that is required once a model is chosen is to calculate the derivatives of the 
potential and no further calculation is necessary. 
It will be useful in the following to define 



Xr = a 2 



/ 8ttG 



(3.2) 



We finally get the Klein-Gordon equation for multiple fields in the flat gauge, from Eq. (2.17), and using all the field 
equations given in the previous section 



Sf 2 i" + 2H5(p 2 / - W 2 S(p 2I + a 2 
16ttG 



K 



8ttG If i tt if 8ttG 

<-' ,<Pk¥i 



Sip- 



2 K 



+ 



+ 



n 



ViVk 



K 



K 



K 



K 



a2U ;Vi V0K 6l PlK + <Poi (^'fK + X k) SifiK 



K 



K 



+ « 2 E 

K,L 



U., 



+ 



8ttG , 



H 



Sip 1K Sip 1L + F (SipiK, 5ip 1K ) = , 



(3.3) 



2 Following on the same lines as detailed in K. A. Malik, D. Seery and K. N. Ananda, arXiv:0712.1787 [astro-ph]. 
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where F [SipiK , S^Pik) contains gradients and inverse gradients quadratic in the field fluctuations and is defined as 
~^ff) Vu,;V~ 2 (X K 8ip 1K + <p' 0K 5<p 1K ') ' ^-V 2 6<pu ^ <p' 0K 6<p 1K 

y~] yo K ^yiK,imV~ 2 51 (^a<5^ia + <Pok s <Pik') * m - Y ( x k S( Pik + <Pok s <Pik') V 2 ^2(p' 0K S<p 1K 

A A 

Voi Y $ l PiK,i8<p 1Kt l + 4A/V~ 2 ^ (Sipix'^StpiK + 5<pik[iS<pik ^ 

A 



+2^ f^\ 2 
H \ H J 



+ 



4ttG 



+ 



4ttG 



1qj 



A' 



A' 



A 



A 



^-V~ 2 jsTrG^ ((5^iajV 2 ^ia-/ + V 2 5<pikV 2 S<pik + <^ia'V 2 <^ia' + ^yk^yk. ') 



2V 2 ^ (Xa-^ik + ^qa^iaO/j X K 5tp 1K + ^ ^oa^ia/ VokSvikj 



A 



A 



A 



A 



(3.4) 



This is the Klein-Gordon equation at second order for multiple scalar fields in terms of the field fluctuations 5<f2i and 
dipu, at second and first order, respectively, in the flat gauge. No slow roll is assumed. 

Note that terms containing inverse gradients and inverse Laplacians never appear on their own. They are always 
either multiplied by terms containing gradients and Laplacians or contain these themselves. Moreover, the order of 
the "positive" gradients is either the same or higher as that of the "negative" or inverse gradients. Hence the large 
scale limit stays well defined and no infrared cut-off has to be introduced by hand. 

As an intermediate step in the derivation of the Klein-Gordon equation above, the equation is given in Appendix 
B in terms of the field fluctuations as in Eq. (3.3), but also with metric potentials 4>i and B\, which allows for a 
more compact presentation. Using Eqs. (2.20) and (2.21) for <fri and B\, respectively, the intermediate Klein-Gordon 
equation, Eq. (Bl), can readily be solved numerically either in real space or in Fourier space (as outlined in the 
following section). However, for analytical studies Eq. (3.3) above seems preferable. 



IV. THE KLEIN GORDON EQUATION IN CLOSED FORM IN FOURIER SPACE 



In order to solve the Klein-Gordon equation it has become standard procedure to work in Fourier space instead of 
real space (see e.g. Ref. [1]). However, at second order things are slightly more complicated than at first order, where 
the shift to Fourier components can be nearly done implicitly, and we shall look at these complications below. 

A perturbation dip is related to its Fourier component 5<pu as 

5ip(n,x l )= J d 3 k5ip k exp (ikiX 1 ) , (4.1) 

where 6tpk = 6(pk(r),k l ), and k l is the comoving wavenumber. We follow the notation of Ref. [1] and use in the 
following the same symbol for the real space and the Fourier space variable, discriminating the two by the respective 
argument, i.e. 5ip(k l ) instead of Sip^, whenever convenient. 



At first order the procedure to transform the Klein-Gordon equation is straight forward: the only change to the 
equation is to replace Sipu by Sipu(k l ) in Eq. (3.1) and the Laplacian by 

V 2 ^-k 2 , (4.2) 

where k 2 = k % ki. At second order things are slightly more complicated, as already pointed out above. Second order 
variables proper, like 5ip2, follow the same rules as at first order and can simply be replaced. However the terms 
quadratic in the first order variables, like Sipi 2 , have now to be treated using the convolution theorem, which is given 
for two functions / and g, both functions of the coordinates x^, by 



[/<4 = / d 3 A/(f>W), 



(4.3) 
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where we introduced d 3 \ = i ^0yrS 3 (k l — p % — q % ). 

In addition to the Laplace-operator we now also have gradients and inverse gradient operators, which in Fourier 
space translate into 



and similarly the inverse Laplacian 



di > iki , d i > i ^2 ) 



-k~ 



(4.4) 



(4.5) 



We have now assembled all the necessary tools to rewrite the Klein-Gordon equation at second order in terms of 
its Fourier components. Rewriting the first part which is independent of gradients is simple and only necessitates in 
Eq. (3.3) the replacement of 6ip2K by its Fourier component, the usage of Eq. (4.2) to replace the Laplacian, similar 
to the first order case, and the application of Eq. (4.3) to terms quadratic in the first order field fluctuation. 

Rewriting the second part involving gradient terms, Eq. (3.4), is more involved, requiring the application of Eq. (4.3) 
and the use of Eq. (4.4) as well as Eqs. (4.2) and (4.5). We leave the "translation" of individual terms in Eq. (3.4) to 
Appendix D, and give here immediately the result, 

v^^i/^E^^rf) + (^) 2 ^ (w - p ^§^) E^^rf)E^<^^) 



Xj_ {4irG\ 2 Pl q l p m q m +p 2 q 2 
H\H) k 2 q 2 



E PoaMk-O**) E (XkSvikW) + (fi'oK^iK'iq 1 )) 



K 



K 



4ttG 



2 / 

K K 



k 2 



1 n J n 



piqpmq 

p 2 q 2 



E (XxStpiictf) + PokMxV)) E (X K S VlK (q l ) + v' 0K 5 VlK '{q 1 )) 



K 



K 



y j 
n 



87rG E ( PI<1 \1 P2 q 2 ^lK{p l )hlK{q l ) - ^ \ Piq Svik'^K-uSW) 



K 



k 2 



(a-kgY tiki Lp 1 p 3 



k 2 



E {Xk5^k{p 1 ) + ¥>ok<W V)) E XKfyuctf) 

P K K / 



(4.6) 



V. SLOW ROLL 

Substantial simplification can be achieved by employing the slow-roll approximation. In particular it will allow us 
to rewrite the source terms of the second order Klein-Gordon equation solely in terms of the field fluctuations Sipu 
themselves, whereas in Eq. (3.3) these terms also includes time derivatives of the first order field fluctuations. 



A. Background and first order 

The slow roll approximation in the background is given by 

<p'oi - Woi - o , E « u ° ■ 



K 



leading to 



8ttG 



3W#,j + a% Vl = , H z = — a z U . 



(5.1) 



(5.2) 
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Making use of Eqs. (5.2) we find that Eq. (3.2) simplifies to 

Xi = 0. 

We introduce the slow-roll parameters [28] 



H 



and 



m 2 



At first order we get the Klein-Gordon equation from Eq. (3.1), imposing slow-roll in the background, 

Sipu" + 2H5<p u ' - V 2 ^!/ + J2 Mi K 5<Pik = o , 



K 



where we defined the mass matrix 

Mir = a 2 U^ lVK - 24wG<p' 0I <p' K . 
Note that r\u is only small in the adiabatic direction, and hence we keep terms including r\u below. 



(5.3) 



(5.4) 



(5.5) 



(5.6) 



(5.7) 



B. Second order 



In the following we shall work to order e|, where the slow parameter ej is defined in Eq. (5.4). We then get the 
Klein-Gordon equation at second order using slow roll and working in Fourier space from Eqs. (3.3) and (4.6), 

5ip 2I "{k i ) + 2HS<p 2I '(k i ) + k 2 S<f2i(k l ) + ]T (a 2 U m - 24irG<p' 0I tp' 0K ) S^ 2K (k l ) (5.6 

K 

,9i<pk<pl + -^-<PbiU,(p KVL j SipiK{p l )S(piL(q l ) H — ^— a 2 (p' QK SipiK(p l ) U, VKV > I S(piK (g 1 ) \ 



8ttG /",,,! 8irG P[ q l 



H 



d 6 X< 



H q 



K 



K 



K 



k 2 



As expected, the simplification of the second order Klein-Gordon equation in the slow roll limit, Eq. (5.8) above, 
compared to the non-approximate versions given in previous sections, Eqs. (3.3) and (4.6), is considerable. 



VI. DISCUSSION AND CONCLUSION 

We have derived the Klein-Gordon equation at second order in the multi-field case in closed form in terms of the 
field fluctuations in the flat gauge. We have given the full equation in real and Fourier space, and also have given a 
simplified version where we used slow-roll. Having calculated the Klein-Gordon equation at second order in closed 
form for multiple fields in Section III allowed us to show that the Klein-Gordon equation also at this order does not 
contain terms that diverge in the large scale limit, irrespective of imposing the slow- roll approximation. Preliminary 
calculations indicate, that this is the case as well if we include vector perturbations. 

If we want to solve the second order Klein-Gordon equation, we have to proceed iteratively: first solving the 
background equations for the fields ifoi and the Hubble parameter H, Eqs. (2.15) and (Eq. (2.18)). Using these 
solutions we can then solve the first order Klein-Gordon equation (3.1) and substitute the results finally into the 
second order equation, (3.3). This procedure is already implemented implicitly in the slow- roll case discussed in 
Section V, using the suitable slow-roll expressions. However, if we are interested in numerical solutions, this iterative 
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scheme has to be explicitly taken into account. The numerical implementation is however straight forward, and despite 
its complexity, second order perturbation theory was shown to be faster numerically than the ATV-formalism [35] . 

However interesting or even exciting we may find the second order Klein-Gordon equation to be, our aim is more 
likely to calculate what observable consequences the evolution of the fields will have and what it can tell us about 
our early universe model, and not the scalar fields themselves. However, once the field dynamics is known, it is easy 
use expressions relating the field fluctuations to observable quantities. A particularly popular such quantity is the 
curvature perturbation on uniform density hypersurfaces, C, at first and second order. How in particular (2 is related 
to the field fluctuations on flat slices, 6(fu and 6<f2i, used in this article, is given in Ref. [29], and its relation to the 
non-linearity parameter / n i is detailed in Ref. [35]. 
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APPENDIX A: USEFUL BACKGROUND RELATIONS 

The Friedmann and the i — j equations in the background, Eqs. (2.18) and (2.19), can be rewritten to give 

*ffl'-T-* a X«''- (A1 » 

Another useful combination is 

W + 2H 2 = 8irGa 2 U . (A2) 



APPENDIX B: KLEIN-GORDON EQUATION AT SECOND ORDER 



It might prove useful to have an intermediate version of the Klein-Gordon equation at second order, that still has 
the metric potentials at first order, whereas the second order metric potentials and the time derivatives of the first 
order variables have been replaced using the field equations given in Section II B. 

U tipiipK 5tp2K + a 2 U, ViVkVl 5(PikS<Pil + 4</>i ^ a 2 U. VlVK 5tpiK — V 2 Sipu J 

K K,L \ K J 



+ 



8irG 



K 



16ttG 



4Vu,fc-Bi,* + *a 2 U m <l>\ + —^-5ip u 'J2 x kS<Pik 



K 



+ 2^ V- 2 |v 2 BiV 2 0i ~ 0i,fciBi, W + SttGJ] (<WV 2 <W + Vik,;<W/) | 



n 



2 [AitG 



K 



K 



^ V 2 J 8ttG^ (s<p 1K ,iV 2 5<p 1K> 1 + V 2 6<p 1K V 2 6<p 1K + 5cp 1K 'V 2 6<p 1K ' + Viir,/<Vi 



K. 



K 

8ttG 



(Bl) 



APPENDIX C: TENSOR COMPONENTS 



1. The metric tensor 
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The metric tensor up to second order for scalar perturbations is 

-a 2 (l + 20i+(/) 2 ) , 



5oo 



1 



g 0l = a z [ Bx + -B 2 



9ii = a 2 [(1 - 20i - 2 ) <% + 2E hij + E 2 , 



and its contravariant form is 



,00 



; [1 - 20! - 2 + 40? _ B Xik B£] , 
Bl + X -Bl - 2B ltk E£ + 2 tyi - 00 Bl 

(1 + 2^i + 02 + 40?) J*i - + £ 2 « - AE^E^l + 80!^^' + B^f) 

Note, Eqs. (CI) and (C4) are without gauge restrictions, i.e. no gauge has been specified. 



(CI) 
(C2) 

(C3) 

(C4) 
(C5) 

(C6) 



2. Energy-momentum tensor in the flat gauge 

The energy-momentum tensor for N scalar fields with potential U(<pi) is then split into background, first, and 
second order perturbations, using Eq. (2.4), as 



1 v — 1 (Q) V + 0± {\)v + 2 ( 2 )" : 

and we get for the components, from Eq. (2.3), at zeroth order 



I (0)0 



K 



K 



at first order 



(lb" " a 2 



and at second order in the perturbations 



5T (°i)o = — \ Y {^'ok 5{ Pik' - v'aK 2< t>i) - 5U\ , 

K 

' K 



(C7) 

(C8) 

(C9) 
(CIO) 

(Cll) 



<5T ( ° 2)0 = — V \<p' 0K 5<P2k' - 4<Pok<I>i s <Pik' - <p' K 2< h + ^K<i>\ + 8<Pik' 2 + a 2 8U 2 

a K L 

+5<p 1K ,i5(p 1Ki l - ip' 0K B 1M B x k , 
ST° 2 )i = -^2 Y WokS<P2k ,i - ±4>w' QK $yiK ,i + 25<p 1K '6<p 1Kii ) , 

K 

ST( 2 )j = ^2 Y [<Pok s <P2k - 4(^o K 0i<5<^ik' - <p' 0K 2 <h + 4:foK 2 4>i + SfiK 2 - S(fi 1K .iS(p 1K l - tp' 0K 2 B ltk B^ 



(C12) 
(C13) 



K 



1-2 

-2ip' QK 8ip 1K ,iB h l - a 2 SU 2 5^ + — (ip' 0K B\ + 6tpi K ) * 6<pi K ,j ■ 



(C14) 
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APPENDIX D: RELATIONS IN FOURIER SPACE 

Here we give the Fourier transform of some of the terms in Eq. (3.4) (integrating out the p l dependence and omitting 
the (27r)~ 3 factor: 

[W] fe = j d 3 q SipitfWMV ~ Q*) , 

[6<pi,i6<p li l ] k = - [d 3 qq i (k i -q i )6 Vl (q i )5^ 1 (k i -q i ), 
[5^,iV- 2 6 Vh l ] k = J d 3 q q- 2 q t (tf - q l )5 Vl {q l )5^(k % - q l ) , 

[V 2 (6<p 1 , ij dr 1 dT 1 6<p 1 )] k = -k 2 Jtfqik-q^qWiki-qJikj-^SMWMV-tf)- ( D1 ) 

The terms of Eq. (3.4) not listed above are similar to the ones given in Eq. (Dl) and can be readily derived using the 
expressions above. 

APPENDIX E: CHANGES COMPARED TO EARLIER VERSIONS OF THE PAPER 

Here we give a list of changes in the second order governing equations compared with the version of the paper 
originally published on the archive, astro-ph/0610864vl. This should make it easier to verify that one is using the 
(most) correct version of the equations . . . However, we do not list mere "cosmetic" changes or rearrangements in the 
equations. 

Changes in version 2: 

1. Eq. (2.24): wrongly set brackets corrected 

2. Eq. (2.29): missing a 2 added 

3. Eq. (3.3): ¥?oi! 12k ^Pik' 2 m 2nd line removed, and rewritten 

4. Imposing slow roll on first order Klein-Gordon equation removed in Section VA (only valid on large scales) 

5. Eq. (5.8): second order Klein-Gordon equation equation in slow roll limit completely rewritten in light of above 
changes 

6. Eq. (Bl): ^^P^k' 2 m square bracket in 2nd line removed 

Changes in version 3 (corresponding to the published version): 

Relaxing the requirement that the slow-roll parameter rju is small (only needs to be small in the adiabatic 
direction) leads to additional terms (of the form tpoU^SipiSifi) in the slow roll Klein Gordon equation Eq. (5.8). 

Changes in version 4 '■ 

Second order field equations corrected (Eqs. (2.26), leading to changes in the second order Klein-Gordon equation, 
Eqs. (Bl), (3.3), (4.6), and (5.8). 

Changes in version 5: 

Second order field equations, Eqs. (2.27), and (2.28) corrected (the concept of an "off-trace" equation proved to be 
not useful at second order), leading to changes in the second order Klein-Gordon equation, Eqs. (Bl), (3.3), (4.6), and 
(5.8). Background slow roll equations in conformal time, Eq. (5.1), corrected leading to changes throughout Section 
V. Also changes in definition of second order flat field perturbation, though leading to no changes in this paper. 
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